United States Patent uq 

Erdol et al. 



US005 103427 A 

[ii] Patent Number: 
[45] Date of Patent: 



5,103,427 
Apr. 7, 1992 



[54] METHOD AND APPARATUS GENERATING 
HIGH RESOLUTION DATA AND ECHO 
IDENTIFICATION 

[75] Inventors: Nurgun Erdol, Highland Beach, Fla.; 

Louis E. Roemer, Akron, Ohio 

[73] Assignee: The University of Akron, Akron, 
Ohio 

[21] Appl. No.: 527,197 

[22] Filed: May 22, 1990 

[51] Int. CI. 5 G03B 42/06 

[52] U.S.C1 367/7; 367/11; 

128/660.01 

[58] Field of Search 367/7, 11, 87, 135; 

364/421; 342/192; 128/660.01; 73/602, 607, 

620, 627 

[56] References Cited 

U.S. PATENT DOCUMENTS 

4,862,180 -8/1989 Martin 342/192 

Primary Examiner— Daniel T. Pihulic 

Attorney, Agent, or Firm — Oldham & Oldham Company 

[57] ABSTRACT 

The present invention is generally directed to an appa- 



ratus and method of echo identification wherein the 
sources of echoes generated from interface boundaries 
in a medium may be more effectively identified in the 
processing of obtained waveform data to generate 
higher resolution in the obtained data. The method of 
the present invention is particularly useful for perform- 
ing ultrasound reflectometry and medical sonography 
by obtaining data representative of a medium being 
studied and suppressing signals from overlapping echo 
waveforms generated from obstacles adjacent an object 
under test. More particularly, the method of the inven- 
tion enables the sources of such echo waveforms to be 
more effectively identified to yield an indication of the 
contribution of overlapping echoes to generate higher 
resolution data signals. The processing includes gener- 
ating a power spectral estimate of the data which is 
utilized to obtain a phase estimate thereof. The Maxi- 
mum Entropy Estimation Method is used to obtain a 
reliable phase estimate which can then be processed to 
yield an indication of the delay time to an interface 
boundary in the medium. The sources of the echo wave- 
forms can then be located and used to increase the reso- 
lution of the obtained data. 



21 Claims, 10 Drawing Sheets 
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METHOD AND APPARATUS GENERATING HIGH 
RESOLUTION DATA AND ECHO 
IDENTIFICATION 

5 

TECHNICAL FIELD 

The present invention is generally directed, to on 
apparatus and method of echo identification wherein 
the sources of echoes generated from interface bound- 
aries in a medium may be more effectively identified in 
the processing of obtained waveform data to generate 
higher resolution in the obtained data. The method of 
the present invention is particularly useful for perform- 
ing ultrasound reflectometry and medical sonography 
by obtaining data representative of a medium being 15 
studied and suppressing signals from overlapping echo 
waveforms generated from obstacles adjacent on object 
under test. More particularly, the method of the inven- 
tion enables the sources of such echo waveforms to be 
more effectively identified to yield an indication of the 20 
contribution of overlapping echoes to generate higher 
resolution data signals. 

BACKGROUND OF THE INVENTION 

In various disciplines such as ultrasound reflectome- 25 
try, medical imaging using sonography, radar, seismol- 
ogy, echocardiography and other similar areas, a com- 
mon problem exists in that waveform data obtained in 
these areas will often include reflections from adjacent 
obstacles which will often overlap the desired echo or 30 
reflected signals. The overlapping echo waveforms 
tend to obscure the desired signals and result in lower 
resolution which inhibits effective and accurate inter- 
pretation of the data. I ultrasound reflectometry which 
is used for medical imaging as well as non-destructive 35 
testing applications, such obscuring echoes result in the 
inability to form high resolution images thereby degrad- 
ing the diagnostic capabilities achievable by the proce- 
dure. Medical imaging has become an important aspect 
of present day medical technology and presently in- 40 
eludes techniques such as X-ray tomography, magnetic 
resonance imaging (MRI) or CAT-scan techniques 
which themselves require relatively high cost, elaborate 
equipment and processing capabilities which inhibit 
their effective use in many situations. If higher resblu- 45 
tion and improved ultrasound images could be obtained, 
this technique may considerably reduce the cost as 
compared to other imaging techniques while providing 
more accurate diagnostics in medicine. 

Ultrasound has been known to be applicable for medi- 50 
cal imaging as well as non-destructive testing as it is a 
mechanical wave phenomena wherein a medium under 
test will enable generated ultrasonic pulses or waves to 
propagate therein. In medical imaging techniques such 
a medium will be the soft tissue of the body with reflect- 55 
ing objects being the internal organs thereof. The wave- 
forms used in ultrasound reflectometry as well as medi- 
cal sonography are typically several cycles in length 
and are rapidly damped to a few wavelengths, at the 
resonant frequency of the transducer. The transducer is 60 
typically a piezoelectric crystal which gives rise to a 
sonic wave propagating into the medium under study 
and reflecting as echoes from reflecting surfaces within 
the medium which are detected by the reverse of the 
piezoelectric effect. In medical sonography for exam- 65 
pie, the transducer normally serves as both a transmitter 
and a receiver, and therefore induced vibrations must be 
damped quickly to avoid a very long ring-down time in 



order to receive echo signals. When low frequencies are 
chosen to obtain low attenuation of the waveforms 
through the medium under test, then the time wave- 
forms of reflections , from adjacent obstacles will often 
overlap. The overlap of echoes makes it difficult to 
interpret where the sources of the echoes are located 
and results in extremely difficult signal processing prob- 
lems to obtain high resolution data which avoids the 
contribution of unwanted overlapping echo informa- 
tion. 

One processing technique being echo envelope pro- 
cessing has been utilized, but it has been found that the 
received information is still significantly contaminated 
by obscuring overlapping echoes. The overlapping 
echo waveform phenomena resulting from the low 
frequencies chosen for ultrasound reflectometry as 
based upon attenuation and resolution constraints have 
been noted in literature such as found in an article by R. 
C. Kemerait and D. G. Childers entitled "Signal Detec- 
tion and Extraction By Cepstrum Techniques" IEEE 
Transactions on Information Theory, Volume 11-18, 
No. 6, pp. 745-749, November, 1972, as well as a publi- 
cation by J. Blitz, entitled "Ultrasonics: Methods and 
Applications", Van Nostrand Reinhold Company, New 
York, 1971. In the first of these publications, there is set 
forth a technique for decomposing a composite signal of 
unknown multiple wavelets which overlap in time. 
Several prior procedures for achieving the decomposi- 
tion of superimposed signals include inverse filtering 
wherein a signal is transformed by a linear time-invari- 
ant system, whose Fourier transform is the reciprocal of 
the transform of the signal components to be removed. 
In such a method, the signal must be known and the 
signal to noise ratio must be quite large. Decision The- 
ory has also been used to decompose superimposed 
signals to estimate the echo amplitude and arrival times, 
but only if the signal wave shape is known. If the wave- 
let waveshape and number of echoes are unknown, 
other techniques have been looked to. The Cepstrum 
techniques of echo detection accomplish decomposition 
by means of a function of the power spectrum of the 
received signals to determine the timing and relative 
amplitudes of echoes in the system. The signal wave- 
form is then extracted by means of complex Cepstrum 
techniques. 

It has also been found that phase information can be 
used in the estimation of time delays between two re- 
ceived signals as described in an article by A. G. Piersol, 
entitled "Time Delay Estimation Using Phase Data", 
IEEE Transactions on Acoustics, Speech and Signal 
Processing, Volume ASSP-29, No. 3, pp. 471-477, June 
1981. Such an estimation of time delays between two 
received signals using phase measurements relied upon 
the use of straightforward regression analysis proce- 
dures on phase estimates at properly selected frequen- 
cies in the frequency domain. Such analysis of the phase 
information was found to yield time delay estimates 
having realistic error assessments based upon non-para- 
metric variance calculations. 

Similarly, a publication by T. F. Quatieri, Jr. and A. 
B. Oppenheim, entitled "Iterative Techniques for a 
Minimum Phase Signal Reconstruction from Phase or 
Magnitude*', IEEE Transactions on Acoustics, Speech 
and Signal Processing, Volume ASSP-29, No. 6, pp. 
1187-1193, December, 1981, develops iterative algo- 
rithms for reconstructing a signal from a partial specifi- 
cation thereof in the time or frequency domains. The 
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technique utilizes iterative algorithms for reconstruct- 
ing a minimum or maximum phase signal from the phase 
or magnitude of its Fourier transform. Additionally, a 
phase unwrapping algorithm is proposed which is im- 
plemented by applying the Hilbert transform to the 5 
logarithmic value of the Fourier transform of a mini- 
mum phase sequence. 

Although the information content of the phase com- 
ponent of systems has been recognized in different ap- 
plications, the phase unwrapping techniques have inhib- 10 
ited their effective use with respect to processing in 
ultrasound reflectometry or similar systems where ob- 
scuring overlapping echoes degrade resolution. For 
example, in a publication by E. Poggiagliomi, A. J. 
Berkhout, M. M. Boone, entitled "Phase Unwrapping, 15 
Possibilities and Limitations", Geophysical Prospect- 
ing, No. 30, pp. 281-291, 1982, a phase unwrapping 
technique is set forth wherein the phase spectrum can 
only be correctly unwrapped between notches in the 
amplitude spectrum and other experimental difficulties 20 
are also present in the described method. 

Some researchers have turned to the principle of 
Maximum Entropy, which involves autoregressive 
modeling, and which has been applied successfully to 
generate power spectral estimations. In publication by 25 
J. P. Burg, entitled "Maximum Entropy Spectral Anal- 
ysis", Modern Spectrum Analysis, edited by D. G. 
Childers, IEEE Press, New York, pp. 34-41, 1978. 
Comparison of the maximum entropy method to some 
traditional techniques was analyzed in a book by A. V. 30 
Oppenheim and R. W. Schafer, entitled "Digital Signal 
Processing, Printice-Hall, Englewood Cliffs, New Jer- 
sey, 1975, wherein the maximum entropy method was 
found to produce higher frequency resolution and yield 
data of better dependability as it depends only upon the 35 
available data and requires simplified storage require- 
ments owing to the infinite impulse response structure 
of the resulting all pole filter. The all pole Filter is found 
to yield a smoothly changing phase estimate, which 
does not appear to have the experimental difficulties 40 
reported with various other methods to obtain phase 
estimates. Such a method has been studied by N. Erdol 
in a PhD dissertation entitled "Use of the Maximum 
Entropy Method for Phase Estimation" at the Univer- 
sity of Akron, Akron, Ohio. In this study, it was found 45 
that phase information obtained by using maximum 
entropy spectral estimates could be utilized to subse- 
quently obtain delay information from the maximum 
entropy phase estimation by means of discrete Fourier 
transformation techniques. 50 

SUMMARY OF THE INVENTION 

Based upon the foregoing, it has been found that 
there is a need to provide an apparatus and method 
which can be used to identify the source of echo signals 55 
to enable higher resolution data to be obtained by mini- 
mizing the contribution of overlapping echoes which 
may obscure the desired echo signals. The present in- 
vention is directed to providing a method and apparatus 
for identifying the locations of the sources of echoes 60 
contributing to data received in systems in applications 
such as ultrasound reflectometry, medical sonography, 
radar/sonar applications, acoustic non-destructive eval- 
uation and testing, seismology, echocardiography and 
the like. The present invention is also directed to a 65 
method and apparatus for improving the resolution of 
received linear frequency modulated signals found in 
such applications^ and particularly to a method and 



apparatus for performing ultrasound reflectometry and 
identifying the source of echoes in this application. 

Generally, the method comprises the steps of obtain- 
ing echo signals in the form of analog data and thereaf- 
ter transforming the analog data signals into a digital 
data sequence. A phase estimate of the digital data se- 
quence is obtained by generating a power spectral esti- 
mate of the data signals by means of a technique known 
as the Maximum Entropy estimation method wherein 
the power spectral estimate is processed to yield the 
phase information. Using phase information, an estimate 
of the delay time of received echo signals which are 
reflected from an interface boundary in a medium can 
be obtained. With an estimate of the delay time, the 
distance to the interface boundary may be determined 
to locate the sources of echo signals, wherein such in- 
formation may then be utilized to minimize the contri- 
bution of overlapping echoes to enhance the resolution 
of the desired echo signals in the system. 

It is therefore a main object of the invention to pro- 
vide a method for echo identification to improve the 
resolution of received echo signals. 

It is another object of the invention to provide for 
high resolution of small structures while retaining low 
attenuation of current signal sources. 

Another object of the invention is to provide an appa- 
ratus and method whereby phase information may be 
obtained from a digital data sequence of reflected echo 
signals to yield an estimate of delay time which can be 
far less than the reverberation time of obscuring echoes 
to minimize reverberation artifacts. 

Yet another object of the invention is to provide a 
method for improving the resolution of linear frequency 
modulated signals by using the method of Maximum 
Entropy at all computational stages which relies only 
upon available data and provides higher frequency reso- 
lution as well as simple storage requirements owing to 
the infinite impulse response structure of the resulting 
all pole filter utilized in the method. 

BRIEF DESCRIPTION OF THE DRAWINGS 

These and other objects of the invention will become 
more readily apparent from a reading of the detailed 
description of the invention in conjunction with the 
drawings, wherein: 

FIG. 1 shows an apparatus for generating and receiv- 
ing ultrasound through an experimental medium, typi- 
fying the apparatus of the invention as used in ultra- 
sound reflectometry; 

FIG. 2 sets forth a simplified representation of an 
ultrasonic pulse signal path within a test medium indi- 
cating the sources of echoes resulting from interface 
boundaries within the medium or obstacle, which may 
result in overlapping echo information; 

FIG. 3 sets forth a block diagram representation of 
the preferred method for echo identification and im- 
proving resolution as associated with a technique such 
as ultrasound reflectometry; 

FIG. 4 sets forth a block diagram setting forth in 
more detail the steps of the method in the present inven- 
tion to identify the sources of echo information; 

FIGS. 5-7 set forth graphical representations of ex- 
perimental results showing a received echo signal with- 
out multiple reflections wherein the amplitude will 
provide digitized values for subsequent processing, the 
phase information obtained therefrom by the method of 
the invention and the estimate of the time delay of the 
received signals respectively; and 
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FIGS. 8-10 set forth graphical representations of 
received echo signals with multiple reflections present 
yielding phase information by means of the method in 
the present invention and an estimate of the time delay 
respectively giving an indication of the source of echo 5 
information to enable estimation of echo contributions 
and to improve resolution of the received signals. 

DETAILED DESCRIPTION OF THE 

INVENTION 10 

Turning now to FIG. 1, an apparatus usable to obtain 
data as weil as to perform processing on the received 
data is shown for use in an ultrasound reflectometry 
method. The invention will be described in terms of 
ultrasound reflectometry but it should be understood 15 
that the apparatus and method may be modified to suit 
other applications such as radar, seismology, echocardi- 
ography and other similar technologies. The apparatus 
as shown in FIG. 1, shows an experimental arrangement 
to show the method of the invention for a simplified 20 
arrangement comprising a water tank 10 having a sim- 
ple reflecting structure 12 therein. The experimental 
arrangement including tank 10 includes a first ultrasonic 
pulse producing transducer 14 arranged on one side of 
the reflecting structure 12 and a second ultrasonic trans- 25 
ducer 16 arranged on the opposite side of the reflecting 
structure 12 to receive echo waveforms. It should be 
recognized that depending upon the particular applica- 
tion of the invention, one or more transducers may be 



an ultrasonic beam is incident upon a boundary between 
two interfaces, part of the beam will be reflected from 
the boundary as shown at 30 while part of the beam is 
transmitted into the second medium shown at. 32 being 
the obstacle of FIG. 2. Within the obstacle, a second 
interface is encountered where again part of the inci- 
dent beam on the second boundary will be reflected as 
shown at 34 while part will be transmitted through the 
boundary into the medium lying beyond the boundary 
as shown at 36. This process continues within the obsta- 
cle creating reverberations wherein the ultrasonic beam 
will be attenuated by absorption and scattering until it is 
ultimately dissipated in the form of random molecular 
motion or heat. As shown in FIG. 2, an incident wave 
upon an obstacle having two boundary surfaces will 
create a number of reverberation echo signals wherein 
each of the echoes will have a different delay time cor- 
responding to the thickness of the obstacle as shown in 
FIG. 2. The fraction of the incident beam intensity 
which is reflected depends upon the acoustic imped- 
ances of the mediums about the boundary, and thus the 
particular medium being tested will directly influence 
the amount of reverberation occurring within the sys- 
tem and the potential of generating additional echo 
waveforms in the system. It is also noted that the frac- 
tion of the ultrasonic beam that remains after reflection, 
which is transmitted into the second medium, will in 
general deviate from its original direction wherein the 
degree of deviation will depend upon the relative veloc- 



utilized. For example, in the application of medical 30 ity of sound in the two media. The direction of the 



sonography, one. ultrasonic transducer will both gener- 
ate ultrasonic pulses and detect echoes that return from 
reflecting surfaces within the body being studied. Typi- 
cally, an ultrasound transducer comprises a thin piezo- 
electric crystal, which is made up of dipoles that are 35 
firmly bound in a crystalline structure. In ultrasound . 
reflectometry including medical sonography, the piezo- 
electric crystal is made to vibrate at its natural fre- 
quency by a single voltage pulse generated by an output 
pulse circuit 18 having power supply 20 and a triggered 40 
output timing circuit 22. The waveforms convention- 
ally used in ultrasound reflectometry are normally sev- 
eral cycles in length, wherein a typical pulse length is 
between 1 to 10 microseconds. Low frequencies are 
normally chosen to obtain low attenuation through the 45 
medium under test and the output pulse is rapidly 
damped to maintain a short ring-down time between the 
initiation of the ultrasonic wave and the cessation of all 
vibrations in the transducer. The ultrasound waveforms 



transmitted beam can be obtained by the application of 
Snell's law similar to light. 

Although a relatively simple reflecting structure is 
shown in the apparatus of FIG. 1, in medical imaging 
the medium through which the ultrasonic pulses will 
travel will be the soft tissue and the reflecting objects 
will be the internal organs. Due to the low frequencies 
used in ultrasound reflectometry to obtain low attenua- 
tion through the medium under test, the time wave- 
forms of reflections from adjacent obstacles, particu- 
larly in more complex systems, will often overlap. The 
overlap of echoes makes it difficult to interpret where 
the sources of the echoes are located, and thus to ac- 
count for their contributions in the received echo sig- 
nals thereby tending to degrade resolution of the re- 
ceived signals. 

Turning to FIG. 3, there is provided a block diagram 
of the apparatus and method in the invention to identify 
the source of echo information. In general, the appara- 



generated by transducer 14 will be detected by trans- 50 tus and method of the invention include means for de- 



ducer 16 in the form of analog data signals which are 
transformed by a conventional A/D converter to en- 
able display on a digital oscilloscope 24 after which the 
digital data signal may be processed by the method of 
the invention using computer 26. It is again reiterated 
that the particular experimental apparatus shown in 
FIG. 1 may be modified in a known manner by those of 
ordinary skill for use of the invention in the other areas 
of technology as previously described, wherein linear 



tecting echo waveform data to be digitized thereby 
forming a digital data sequence as indicated at 50. The 
digital data signal is thereafter processed at 52 to obtain 
a phase estimate of the data signal containing informa- 
55 tion correlating to the time delay of the echo wave- 
forms defected in the system. At 54, means are provided 
to obtain an estimate of the delay time of the echo wave- 
form data from the phase estimate obtained at 52. Fi- 
nally, at 56 there are means to obtain a distance to a 



frequency modulated signals generated from echo 60 reflecting boundary being the source of the echo wave- 



waveform data may be obtained by suitable detection 
means and subsequently processed to yield an estimate 
of the echo contributions and to identify the source of 
such echoes to yield high resolution data or images. 

In FIG. 2, there is shown the signal path for an inci- 
dent wave through an obstacle in its path giving some 
idea of echo waveform contributions for a simple re- 
flecting structure in the path of an incident wave. When 



65 



form data which correlates to the delay time calculated 
in step 54. These steps form a very generalized descrip- 
tion of the method of the invention as a means of echo 
identification of received waveform data. 

The invention will now be described in more detail 
relative to the general procedures as outlined with ref- 
erence to FIG. 3. Turning to FIG. 4, the method of the 
invention includes a first step in computing the coefflci- 
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ents of a minimum phase filter H(z) to form a prediction 
filter 



(i) 



H{2) 



10 



which adequately represents the digital data signal x(n). 
Prediction filters are computed at 62 for the even and 
odd data sequences generated at 60. The design of the 
prediction filter, represented by the even and odd con- 
stituents of the data sequence x(n), is particularly impor- 
tant in applications such as ultrasound reflectometry 
and other similar environments where the digital data 
signal x(n) is only partially available due to the nature of 
the data and its method of acquisition. In such systems, 15 
the unknown portion of x(n) must be predicted and the 
power spectral density is necessarily estimated, which 
require the predictive filter to yield as accurate a pre- 
diction of the unknown portion of the data as possible. 
As the power spectral density can be computed only if 20 
all of the data x(n) is known, the two methods including 
the prediction of the unknown portion of x(n) and the 
estimation of the power spectral density are related and 
will benefit from a model which uses only the available 
data to accomplish this task. 25 

Modeling the system as an all pole, minimum phase 
filter provides a linear model, but it is noted that the 
processing preformed in the method is not linear. The 
all pole, minimum phase filter also simplifies storage 
requirements owing to the infinite impulse response 30 
structure of the filter. 

The unit sample response y(n) of Y(z) is not directly 
related to the data x(n), however it can be shown that 
their respective autocorrelation values are related by: 



/? ja (n)=^n)n=0,l, . . . , N 



(2) 



wherein N is the order of the all pole, minimum phase 
filter. It thus can be seen that R>y(n), for n greater than 
N, provides a non-zero extension to Rxx(n). Since the 
power spectrum is the Fourier transform of its autocor- 
relation sequence, it can thus be shown that 



35 



40 



(3) 



1 7— c/w 



45 



50 



is an estimate of the power spectrum of the data x(n). 
The estimation of the power spectrum is independent of 
the phase of x(et> w l), and the error noise power of the 
predictive filter is <r 2 . 

The invention utilizes the above mentioned filter 
design methodology to determine the distance to a 
source of echo waveform data such as the distance to a 
tissue interface in a medical sonography application. A 
brief description of the mathematical model to deter- 55 
mine this distance is shown in the following. 

A wide-band, short time duration pulse x(t), and its 
reflections x from an interface boundary at a distance d,* 
can be represented at a receiver as: 

60 

xK/)=*(0+ap<'-'l)+a2*('-'2)+ (4) 

where the delay time of the waveform data from the 
interface boundary, t</=2dyv, being the round trip 
distance to the interface boundary divided by the speed 65 
of the reflections in the medium. The term, x, represents 
the widening of the echo waveform as it travels through 
the medium and a/ accounts for the attenuation of the 



waveform in the medium. As mentioned before, short 
distances involved and the widening of the echo wave- 
form result in echoes which overlap. Windowing tech- 
niques are used to remove the overlapping pans of the 
received waveforms, but are often found to severely 
truncate the reflected waveform. The signal to be 
treated after windowing can be represented as: 



^/-fl-€)=aiJt(/-ri) »rO-/|-e) 



(5) 



where r is the width of window ,w, being symmetrical 
about ti+e, which is expected to be the center of data 
of the reflected signal. Due to the effects of dispersion in 
the medium, the center of the reflected signal is not as 
well defined as the incident signal, and it is thus the goal 
to compute ti=2di/v to enable computation of the 
distance d[. 

It is known that the incident signal x(t) is symmetrical 
and it is reasonable to model the medium through which 
the echo waveform travels as a linear phase system. 
Translating the center of the window to the origin 
(t=0) causes the center of the reflected signal to be e 
away from the origin, where e may be positive or nega- 
tive. Since a translation by e in time implies a linear 
phase term in the frequency domain, the problem be- 
comes one of finding the linear phase component of 
aix(t— e) from its asymmetrically truncated form 
y(t-€). 

The Maximum Entropy spectral estimation philoso- 
phy and methodology is utilized to estimate the phase 
and thus the delay time and therefore the distance to the 
interface boundary being a source of echoes. The Maxi- 
mum Entropy Method enables finding an accurate esti- 
mation of the phase of a signal which would normally 
require complete knowledge of the signal, but wherein 
such knowledge is unavailable due to windowing. Al- 
though the Maximum Entropy Method uses a non-lin- 
ear criterion of goodness, the predictive filter which is 
computed is a linear filter which generates a substitute 
signal for the truncated one caused from windowing. 
This new signal is not truncated but is not a point by 
point equivalent of the truncated one, although provid- 
ing an accurate estimation. The power spectrum of the 
substitute signal has been found to be very close to the 
power spectrum of the original signal before truncation, 
and thus addresses the problem of using a power spec- 
tral estimation method to obtain phase estimation. 

For simplicity of notation, we will use x(t) to denote 
the signal whose phase we wish to estimate. The signal 
x(t) is sampled to give a number of equally spaced sam- 
pled data points x(n) whose Fourier transform is X(ef a ). 
At 66 in FIG. 4, there are provided means to obtain the 
phase estimate of the sample data for the process x from 
the known data samples x(0), x(l), ...» x(M), where M 
is the data length being an even integer. To obtain the 
phase estimate, the time series x(k) is written as the sum 
of an even series, e(k), and an odd series o(k), where the 
point of symmetry is chosen at N5 =M/2. The power 
spectra of the time series and its even and odd compo- 
nents are related by 



(6) 



wherein the Fourier transforms of each may be repre- 
sented by 



0) 
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with the even and odd components having real and ded in quotation marks because it has the unusual di- 

imaginary portions. The phase 6(a)) may be defined as merision of time incorporated therein. 

Once the delay id is known, the distance d to the 

0^ (g) boundary may be determined at 70 which provides 

$(&>) = Arc tan 5 means to estimate echo contributions by determining 

the source of the echo waveform data, 
where 0(e/*>) and E(e/») are the real functions of o for Turning now to the experimental results shown in the 
the even and odd components and the factor e-^is reinaining Figures, the apparatus and method of the 
due to the fact that the origin of symmetry is T^N^Ato invention were applied to a simplified experimental 
where N, is the number of data samples to the center of 10 model to give some indication of the improved resolu- 
the data block and Ato is the sample spacing. From this, tion and method of identifying the location of the 
it may be shown that sources of echoes in the sampled data. In the experimen- 
tal model, a sheet of aluminum having a thickness 0.16 
c jj*>\ c (9) 0111 was placed into the path of ultrasound similar to that 
o»I2tf(«)J = *T " jj 15 shown in FIG. 1. In the apparatus of FIG. 1, the ultra- 
s'©*") + SoA* 0 *) sound transducers 14 and 16 were 3.5 MHz focused 

transducers, mounted on opposing sides of a plastic tank 

wherein the even and odd power spectral terms SMC") containing water. In the analog to digital conversion of 

and S»(e>) may be expressed in terms of the even and the received data were oversampled at 20 M 

odd components of the predictor polynomial H(z) 20 samples/second and then averaged for an effective sam- 

y* eldin S pie rate of 10 M samples/second. Although the trans- 
ducer energy is centered around the 3.5 MHz resonance 

Pptf^ofrWe-/") - PoHt/to&WJfi-**) < 10 > of the transducer, the information content in the signal 

Cos[2e(o>)J = p^^Hjp-j**) + PoHdofrWtte-**) ' phase and amplitude is distributed around the unit circle 

in the z plane. It has been found that the poles of the all 

To be consistent with the principles of Maximum pole predictive filter convey their information as the 

Entropy, the point of even and odd symmetry is chosen unit circle is traversed, wherein the most rapid changes 

at the median of the available data at N 5 =M/2. The occur in the vicinity of the poles. Although the signal 

even series, e{k) and the odd series o(k), are put through 3Q energy will be concentrated near the 3.5 MHz trans- 

their respective Maximum Entropy filters wherein the ducer frequency, the computation of the prediction 

outputs will be the power spectra of ejt, ojt accordingly. error filter coefficients, and therefore the pole locations, 

To obtain the cosine of twice the phase of x*, the ratio will evidence the information from the poles of the 

of the difference and the sum of the even and odd series predictive filter on the whole unit circle. If only the 

power spectra are computed. 35 information on the unit circle adjacent to the poles 

Since the Maximum Entropy filters through which would be to window the result, thereby eliminating 

the even and odd series are processed are all pole, mini- useful information for the Maximum Entropy phase 

mum phase filters, the right hand side of Equation 10 estimation. Therefore, the entire unit circle is traversed 

converges for all values of <o, and the phase 6 (&>) can be the ,Z plane using transform calculations to obtain the 

calculated from Equation 10. ^ entire information content of the waveform data for the 

After the phase estimate of the sample data has been phase estimation using the Maximum Entropy Method, 

obtained by the Maximum Entropy Method, this esti- Further in the experimental configuration, a bi-junc- 

mate is used at 68 which provides a means to calculate $ on transistor was used to discharge a capacitor across 

the time delay of received echo waveforms by obtaining ^ transmitter transducer 14 to obtain a sufficiently 

the power spectrum of the phase estimate, again using 4J large The received signal at transducer 16 was 

the Maximum Entropy Method. As an example, a pulse sufficiently large to be digitized by a Tektronix Model 

of ultrasonic energy traveling through a uniform loss- 2 220 digital oscilloscope at 24 without prior amplifica- 

less medium will be reflected from a boundary and the ^ 

echo will be received at the point of origin of the signal. In a flfSt examplej a direct transmission of the pulse 

The incidence signal, x(t) of ultrasonic energy will from the 3.5 MHz transducer 14 to receiving transducer 

travel a total distance 2d through the medium which has 16 of ^ &ame freque ncy is shown, wherein the obstacle 

a velocity of propagation v. The received signal y(t), u $ removcd ffom ^ water ^ 10 m ^ apparatus of 

will be a function of the delay time, W, where w=2d/v. ^ ^ Tmshi tQ FJG 5 ^ ^ ^ m obsta . 

By preprocessing the echo signal according to the Max- ultrasound path is shown which does not have 

unum Entropy Method to yield the phase estimate de- ^ reflections ^ litude of ^ re _ 

scnbed at 66, it may be shown that ^ wavefonn ^ the digiti2ed value> which My * 

a^a^caxiuw (ID regarded as a relative value. Utilizing the method of the 

invention as previously described, the phase estimate 
where the argument will be proportional to the delay *(<■>) can be calculated from Equation 10 which gives 
time of the echo with all data referenced to the origin of 60 the cosine of twice the phase angle being shown m FIG. 
symmetry being the median of the available data ac- 6. Proceeding with the method as previously defined, 
quired from the reflected echo waveform. The amount the spectrum of the phase estimation data as shown in 
of delay with respect to a known point, being the point FIG. 6 yields the echo identification giving a delay time 
of even and odd symmetry, is captured as the "fre- as shown in FIG. 7. 

quency" of a sinusoid. The "frequency" of the sinusoid 65 Turning to FIG. 8, after the thin aluminum sheet is 
can be extracted by obtaining the Maximum Entropy placed in the water of tank 10 as an obstacle shown at 12 
power spectrum of the phase estimate of the original in FIG. 1, the received signal shows a long, ill-defined 
data or COS28. The term "frequency" has been embed- reflected signal indicating multiple reflections as con- 
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trasted with FIG. 5. The extended echo of FIG. 8 ap- 
pears without clear delineation between the different 
components of the multiple echo, making interpretation 
of the obtained waveform data extremely difficult. 
Using the method of the invention, the phase estimation 5 
being the cosine of twice the phase angle is shown in 
FIG. 9, wherein the oscillatory nature of the curve 
enables identification of the multiple echoes as shown in 
FIG. 10. In FIG. 10, the position of the first peak corre- 
sponds to the time of transit through the aluminum 
sheet which comprises the test obstacle. This transit 
time is the round trip or twice the path length within the 
0.16 cm obstacle. The distance, d, to the boundary may 
be computed from the time delay calculation as shown J5 
in FIG. 10, wherein for the experiment described above, 
this sample calculation is 



</=J(6420m/i) (0.5iu) = 1.6mm 



20 



which is seen to verify and identify this source of multi- 
ple echoes. It is also noted that the relative amplitudes 
of the multiple echoes are consistent. Based upon this 
simplified experimental arrangement, the computed 
value of U was found to be quite accurate and only 2 s 
limited by the resolution obtained in digital transforma- 
tion of the data. The resolution can be improved by 
increasing the number of samples, but at the expense of 
increased filter lengths in the processing. Further exper- 
imental results show that the method provides stable 30 
delay detection with respect to variation of data length 
and additive noise, with accuracy being a function of 
filter length only. Although a simplified model has been 
used to give an indication of results obtainable by the 
method and apparatus of the invention, it should be 35 
apparent that the invention can be extended to more 
complex systems in a variety of areas such as medical 
sonography, ultrasound non-destructive testing, 
echocardiography, underwater sonar, radar and geo- 
physical data as previously described. The invention 40 
may also be applied to dispersive media to detect group 
delays as well as phase delays beneficially. 

Although the invention has been described with re- 
spect to ultrasound reflectometry in some detailed as- 
pects thereof, it should be obvious that many modifica- 
tions and variations of the present invention are possible 
in the light of the general concepts of the disclosed 
invention. It is therefore understood that the invention 
may be practiced other than as specifically described 5Q 
herein and is only limited by the scope of the appended 
claims. 

What is claimed is: 

1. A method for performing waveform processing 
comprising the steps of: 55 
obtaining waveform signals in the form of analog data 

signals by means of detection means, 
transforming said analog data signals to a digital data 

signal, 

obtaining a phase estimate of said digital data signal 60 
by generating a power spectral estimate of said 
digital data signal by means of a technique known 
as the Maximum Entropy Estimation Method, and 
processing said power spectral estimate, to yield 
said phase estimate, 65 

obtaining an estimate of the delay time of said wave- 
form signals which are reflected from an interface 
boundary in a medium, 



detennining the distance to said interface boundary 
by means of said estimate of the delay time to lo- 
cate the sources of said waveform signals, 

processing said waveform signals to account for the 
sources generating these signals to minimize the 
contribution of predetermined waveform signals. 

2.. A method for performing waveform processing as 
in claim 1, wherein, 

the step of obtaining said phase estimate of said digital 
data signal includes modeling the system as an all 
pole, minimum phase filter which yields a smoothly 
changing phase estimate of said digital data signal. 

3. A method for performing waveform processing as 
in claim 1, wherein, 

said step of obtaining a phase estimate of said digital 
data signal includes the steps of utilizing the even 
and odd components of said digital data signal and 
computing a prediction filter for each of said even 
and odd components from which said power spec- 
tral estimate can be generated. 

4. A method for performing waveform processing as 
in claim 1, wherein, 

said step of obtaining a phase estimate of said digital 
data signal includes computing the coefficients of a 
minimum phase filter to generate a prediction filter 
which adequately represents said digital data signal 
in systems where said digital data signal is only 
partially available and the unknown portion of said 
digital signal must be predicted. 

5. A method of performing waveform processing as 
in claim 1, wherein, 

said step of obtaining a phase estimate includes gener- 
ating a power spectral estimate of the even and odd 
components of said digital data signal wherein the 
phase, 0(0)% is related to the power spectral esti- 
mate as represented by the following: 



Cos[20<oO] ■ 



which may be shown to yield 



45 



Cos[2fl(«)l . 



wherein the phase Q{o) can then be calculated. 

6. A method for performing waveform processing as 
in claim 5, wherein, 

said step of obtaining an estimate of the delay time 
utilizes the phase estimate which is related to the 
delay time as follows: 

Cas20(a>)= CosQxatd) 

where the argument is proportional to the delay 
time of the echo waveform data. 

7, A method for performing waveform processing as 
in claim 6, wherein, 

said step of determining the distance to said interface 
boundary utilizes the estimate of the delay time 
wherein the distances related to the delay time by 
the following: 
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where d is the distance to said interface boundary 
and v is the velocity of the waveform signal in the 
medium through which it travels. 

8. A method for performing waveform processing as 
in claim 1, wherein, 5 

said waveform signals are echo signals and said step 
of determining the distance to said interface bound- 
ary enables identifying the echo location of the 
source of echoes wherein said waveform data in- 
cludes overlapping echoes which tend to obscure 10 
waveform data, wherein identification of the echo 
location of the source of said overlapping echos 
allows minimization of the contribution of said 
overlapping echos. 

9. A method for performing waveform processing as 15 
in claim 1, wherein, 

said step of obtaining an estimate of the delay time 
includes obtaining a power spectrum of the phase 
estimate by means of the Maximum Entropy Esti- 
mation Method, wherein the power spectrum of 
the phase estimate curve exhibits a peak at a time 
corresponding to said delay time. 

10. A method for performing waveform processing as 

in claim 1, wherein, 25 
said waveform data is generated in an ultrasound 
reflectometry system. 

11. A method for performing waveform processing as 
in claim 1, wherein, 

said waveform data is generated in a medical sonog- 3Q 
raphy application. 

12. A method for performing waveform processing as 
in claim 1, wherein, 

said waveform data is generated in an acoustic non- 
destructive evaluation and testing application. 33 

13. A method for performing waveform processing as 
in claim 1, wherein, 

said waveform data is generated in a radar system. 

14. A method for performing waveform processing as 

in claim 1, wherein, 40 
said waveform data is generated in a sonar system. 

15. A method for performing waveform processing as 
in claim 1, wherein, 

said waveform data is generated in a seismology ap- 
plication. 45 

16. A method for performing waveform processing as 
in claim 1, wherein, 

said waveform data is generated in an echocardiogra- 
phy application. 

17. A method for performing ultrasound reflectome- 50 
try comprising the steps of: 

generating ultrasound waveforms within a medium 
by a transducer means, 

detecting ultrasonic echo signals in the form of ana- 
log data signals by transducer means, 55 

transforming said analog data signals to a digital data 
signal, 

obtaining a phase estimate of said digital data signal 
by generating a power spectral estimate' of said 
digital data signal by means of a technique known 60 
as the maximum entropy estimation method, and 
processing said power spectral estimate to yield 
said phase estimate, 



obtaining an estimate of the delay time of said ultra- 
sound echo signals which are reflected from an 
interface boundary within said medium, 

determining the distance to said interface boundary 
by means of said estimate of the delay time to lo- 
cate the sources of said echo signals, 

processing said echo signals to account for the 
sources generating the ultrasound echo signals 
wherein the contribution of predetermined echo 
signals is minimized. 

18. A method for performing ultrasound reflectome- 
try as in claim 17, wherein, 

the echo waveform data is generated in a medical 
sonography technique. 

19. A method for performing ultrasound reflectome- 
try as in claim 17, wherein, 

the echo waveform data is generated in a acoustic . 
non-destructive evaluation and testing application. 

20. A method of identifying the echo location of the 
source of echoes in echo waveform data, which com- 
prises the steps of: 

obtaining echo waveform signals in the form of ana- 
log data signals by means of detection means, 

transforming said analog data signals to a digital data 
signal, 

obtaining a phase estimate of said digital data signal 
by generating a power spectral estimate of said 
digital data signal by means of a technique known 
as the Maximum Entropy Estimation Method, and 
processing said power spectral estimate to yield 
said phase estimate, 

obtaining an estimate of the delay time of said echo 
signals which are reflected from an interface 
boundary in a medium, 

determining the distance to said interface boundary 
by means of said estimate of the delay time to lo- 
cate the sources of said echo signals. 

21. An apparatus to perform waveform processing, 
comprising: 

means to generate waveforms in a medium which will 
be reflected for interface boundaries within said 
medium to form said echo waveforms in said me- 
dium, 

means to detect said echo waveforms to generate 
analog echo data signals, 

means to transform said analog echo data signals to a 
digital data signal, 

means for processing said digital data signal to gener- 
ate a power spectral estimate thereof by means of 
maximum entropy estimation, 

means to calculate a phase estimate of said digital data 
signal using said power spectral estimate, and cal- 
culating an estimate of the delay time of said echo 
signals being the travel time within said medium to 
said interface boundary, 

means to calculate the distance to said interfaced 
boundary using said estimate of the delay time so as 
to locate the echo source of said echo waveforms, 

means for processing said digital data signal to ac- 
count for the sources of said echo waveforms 
wherein the contribution of predetermined echo 
waveforms is minimized. 



65 
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[57] ABSTRACT 

A constrained-envelope digital-communications transmitter 
circuit (22) in which a binary data source (32) provides an 
input signal stream (34), a phase mapper (44) maps the input 
signal stream (34) into a quadrature phase-point signal 
stream (50) having a predetermined number of symbols per 
unit baud interval (64) and defining a phase point (54) in a 
phase-point constellation (46), a pulse-spreading filter (76) 
filters the phase -point signal stream (50) into a filtered signal 
stream (74), a constrained-envelope generator (106) gener- 
ates a constrained-bandwidth error signal stream (108) from 
the filtered signal stream (74), a delay element (138) delays 
the filtered signal stream (74) into a delayed signal stream 
(140) synchronized with the constrained-bandwidth error 
signal stream (108), a complex summing circuit (110) sums 
the delayed signal stream (140) and the constrained- 
bandwidth error signal stream (108) into a constrained- 
envelope signal stream (112), and a substantially linear 
amplifier (146) amplifies the constrained-envelope signal 
stream (112) and transmits it as a radio-frequency broadcast 
signal (26). 

29 Claims, 4 Drawing Sheets 
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CONSTRAINED-ENVELOPE DIGITAL- 
COMMUNICATIONS TRANSMISSION 
SYSTEM AND METHOD THEREFOR 

TECHNICAL FIELD OF THE INVENTION 

The present invention relates generally to the field of 
digital communications. More specifically, the present 
invention relates to the field of constrained-envelope digital 
transmitter circuits. 

BACKGROUND OF THE INVENTION 

A wireless digital communications system should ideally 
refrain from using any portion of the frequency spectrum 
beyond that actually required for communications. Such a 
maximally efficient use of the frequency spectrum would 
allow the greatest number of communications channels per 
given spectrum. In the real-world, however, some spectral 
regrowth (i.e., increase in spectral bandwidth) is inevitable 
due to imperfect signal amplification. 

In wireless communication systems various methodolo- 
gies have been used to minimize spectral regrowth. Some 
conventional methodologies utilize complex digital signal 
processing algorithms to alter a digitally modulated trans- 
mission signal in some manner conducive to minimal spec- 
tral regrowth. Such complex algorithmic methodologies are 
well suited to low-throughput applications, i.e., those less 
than 0.5 Mbps (megabits per second), such as transmission 
of vocoder or other audio data. This is because the low 
throughput rate allows sufficient time between symbols for 
the processor to perform extensive and often repetitive 
calculations to effect the required signal modification. 
Unfortunately, high-throughput applications, i.e., those 
greater than 0.5 Mbps, such as the transmission of high- 
speed video data, cannot use complex processing algorithms 
because the processing power required to process the higher 
data rate is impractical. 

A digital signal processing methodology may be used 
with the transmission of burst signals. With burst 
transmissions, the interstitial time between bursts may be 
used to perform the necessary complex computations based 
upon an entire burst. This methodology is not practical when 
continuous (as opposed to burst) transmission is used. 

A conventional form of post-modulation pulse shaping to 
minimize spectral bandwidth utilizes some form of Nyquist- 
type filtration, such as Nyquist, root-Nyquist, raised cosine- 
rolloff etc. Nyquist-type filters are desirable as they provide 
a nearly ideal spectrally constrained waveform and negli- 
gible inter-symbol interference. This is achieved by spread- 
ing the datum for a single constellation phase point over 
many unit baud intervals in such a manner that the energy 
from any given phase-point datum does not interfere with 
the energy from preceding and following phase-point data at 
the appropriate baud-interval sampling instants. 

The use of Nyquist-type filtration in a transmission circuit 
produces a filtered signal stream containing a pulse wave- 
form with a spectrally constrained waveform. The degree to 
which a Nyquist-type pulse waveform is constrained in 
bandwidth is a function of the excess bandwidth factor, a. 
The smaller the value of a, the more the pulse waveform is 
constrained in spectral regrowth. It is therefore desirable to 
have the value of a as small as possible. However, as the 
value of a is decreased, the ratio of the spectrally con- 
strained waveform magnitude to the spectrally uncon- 
strained waveform magnitude is increased. The spectrally 
unconstrained waveform is the waveform that would result 
if no action were taken to reduce spectral regrowth. Typical 
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designs use a values of 0.15 to 0.5. For an exemplary a value 
of 0.2, the magnitude of the spectrally constrained waveform 
is approximately 1.8 times that of the unconstrained wave- 
form. This means that, for a normalized spectrally uncon- 

5 strained waveform magnitude power of 1.0, the transmitter 
output amplifier must actually be able to provide an output 
power of 3.24 (1.8 2 ) to faithfully transmit the spectrally 
constrained waveform. This poses several problems. 
When the transmitter output amplifier is biased so that the 

10 maximum spectrally unconstrained waveform (1.0 
normalized) is at or near the top of the amplifier's linear 
region, all "overpower" will be clipped as the amplifier 
saturates. Such clipping causes a marked increase in spectral 
regrowth, obviating the use of Nyquist-type filtration. 

15 When the transmitter output amplifier is biased so that the 
maximum spectrally constrained waveform (1.8 
normalized) is at or near the top of the amplifier's linear 
region, the spectrally unconstrained waveform is at only 56 
percent (i.e., Via) of the amplifiers peak linear power. This 

20 results in an inefficient use of the output amplifier. 

Also, the biasing of the transmitter output amplifier so that 
the spectrally constrained waveform is at or near the top of 
the amplifier's linear region requires that the output ampli- 

25 fier be of significantly higher power than that required for 
the transmission of a spectrally unconstrained waveform. 
Such a higher-power amplifier is inherently more costly than 
its lower-power counterparts. 

3Q SUMMARY OF THE INVENTION 

It is an advantage of the present invention that a circuitry 
and a methodology are provided that allow a transmitter 
output amplifier to be biased so that the spectrally uncon- 
strained waveform is at or near the top of the amplifier's 

35 linear region without incurring clipping of a spectrally 
constrained waveform. 

It is another advantage of the present invention that a 
circuitry and methodology are provided that allow a spec- 
trally constrained waveform to have approximately the same 

40 magnitude as the spectrally unconstrained waveform with- 
out effecting a significant increase in spectral regrowth. 

It is another advantage of the present invention that a 
circuitry and methodology are provided which allow a 

45 spectrally constrained waveform to be utilized with a con- 
tinuous transmission scheme. 

It is another advantage of the present invention that a 
circuitry and methodology are provided which allow effi- 
cient use of a transmitter output amplifier, thus allowing 

50 higher power output for a given output amplifier and a given 
bandwidth constraint than was previously feasible. 

It is another advantage of the present invention that a 
circuitry and methodology are provided which allow effi- 
cient use of a transmitter output amplifier, which allows 

55 allowing a lower-power amplifier to be used for achieving 
given bandwidth constraints than was previously feasible, 
thus effecting a significant saving in the cost thereof. 

These and other advantages are realized in one form by a 
constrained-envelope digital communications transmitter 

60 circuit. The transmitter circuit has a pulse -spreading filter 
configured to receive a quadrature phase-point signal stream 
of digitized quadrature phase points and produce a filtered 
signal stream, which filtered signal stream exhibits energy 
corresponding to each phase point spread throughout a 

65 plurality of baud intervals. The transmitter circuit also has a 
constrained-envelope generator coupled to the pulse- 
spreading filter and configured to produce a constrained- 



12/30/04, EAST. version: 2.0.1.4 



6,104,761 

3 4 

bandwidth error signal stream. The transmitter circuit also Digital communications system 20, as depicted in FIG. 1, 
has a combining circuit coupled to the pulse-spreading filter includes a transmitter circuit 22 and a transmitter antenna 24 
and to the constrained-envelope generator, which combining together configured to modulate and transmit a radio- 
circuit is configured to combine the filtered signal stream frequency (RF) broadcast signal 26 to a receiver antenna 28 
and the constrained-bandwidth error signal stream to pro- 5 and a receiver circuit 30, together configured to receive and 
duce a constrained-envelope signal stream. The transmitter demodulate RF broadcast signal 26. Those skilled in the art 
circuit also has a substantially linear amplifier with an input will appreciate that the embodiment of system 20 depicted 
coupled to the combining circuit. is a simplistic one for purposes of discussion only. In normal 
These and other advantages are realized in another form use » system 20 would likely be a complex system consisting 
by a method for the transmission of a constrained -envelope io of many more components and broadcast signals. It will be 
communications signal in a digital communications system. appreciated that the use of such a complex communications 
The transmission method includes the step of filtering a system for system 20 in no way departs from the spirit of the 
quadrature phase-point signal stream to produce a filtered present invention or the scope of the appended claims, 
signal stream, which filtering step spreads energy from each Transmitter circuit 22 has a binary data source 32 pro- 
phase point over a plurality of baud intervals. The transmis- 35 viding a binary input signal stream 34. Binary data source 32 
sion method also includes the step of generating a may be any circuitry, device, or combination thereof pro- 
constrained -bandwidth error signal stream from the filtered ducing input signal stream 34. Input signal stream 34 is 
signal stream and a threshold signal. The transmission made up of binary data that may be pre-encoded in any 
method also includes the step of combining the filtered desired manner. That is, input signal stream 34 may be made 
signal stream and the constrained-bandwidth error signal 20 U p 0 f da t a that has no encoding, concatenated encoding, 
stream to produce a constrained-envelope signal stream. The Reed-Solomon block encoding, or any other form of encod- 
transmission method also includes the step of linearly ampli- ing desired for or required of the communications scheme in 
fying the constrained -envelope signal stream to produce the use. 

constrained-envelope communications signal. The transmis- In the preferred embodiment, input signal stream 34 is a 

sion method also includes the step of transmitting the 25 steam of contimlous data (as contrasted with burst data) 

constrained-envelope communications signal. passing to an input of a convolutional encoder 36 . Convo . 

BRIEF DESCRIPTION OF THE DRAWINGS lutional encoder 36 convolutionally encodes (e.g., Viterbi 

encodes) input signal stream 34 into an encoded signal 

A more complete understanding of the present invention stream 38 ^ ^ of convolut ion a l encoder 36 in trans- 
may be derived by referring to the detailed description and 30 miUer ckcuit 22 &nd a ^ ^0^^ decoder (not 
claims when considered in connection with the Figures, shown) m receiver circuit 30 significantly reduces the error 
wherein like reference numbers refer to similar items rate of lhe overall signal in a manner weU understood by 
throughout the Figures, and: lhose sJdlled m tne art However, convolutional encoder 36 

FIG. 1 depicts a simplified block diagram of a digital may be omitted, 

communications system in accordance with a preferred 35 i Dt erleaver 40 temporally decollates encoded signal 

embodiment of the present invention; stream 38 to produce an mterleaved signal stream 42 ^ 

FIG. 2 depicts a block diagram of a constrained-envelope is, the symbols making up the binary signal stream are 

digital communications transmitter circuit in accordance temporally decorrelated (i.e., separated) in transmitter cir- 

with a preferred embodiment of the present invention; cuit 22 and temporally correlated in receiver circuit 30. This 

FIG. 3 depicts a 16-P-APSK constellation illustrating a is done so that correlated errors produced by downstream 

locus of a quadrature phase -point signal stream over twelve transmitter components, discussed hereinbelow, will then be 

exemplary consecutively mapped phase points in accor- decorrelated through a complimentary de-interleaver located 

dance with a preferred embodiment of the present invention; in receiver circuit 30 before convolutional decoding in 

FIG. 4 depicts a plurality of signal streams in accordance 45 receiver circuit 30. 

with a preferred embodiment of the present invention; In the preferred embodiment, interleaved signal stream 42 

FIG. 5 depicts the phase-point constellation of FIG. 3 passes to an input of a phase mapper 44. Those skilled in the 

illustrating an exemplary locus of a filtered signal stream art wiil appreciate that interleaver 40 is not desired in all 

over the twelve consecutively mapped phase points of FIG. embodiments of transmitter circuit 22, for example when 

3 in accordance with a preferred embodiment of the present 50 convolutional encoder 36 is omitted. When interleaver 40 is 

invention; omitted, encoded signal stream 38 is passed directly to the 

FIG. 6 depicts a pair of Nyquist-type data bursts in input of phase mapper 44. When both convolutional encoder 

accordance with a preferred embodiment of the present 36 and interleaver 40 are omitted, binary input signal stream 

invention; and passes directly to the input of phase mapper 44. 

FIG. 7 depicts a noise-influenced constellation illustrating 55 FIG. 3 depicts a sixteen phase-point polar amplitude and 

constrained-envelope phase-point probabilities of the phase P hase shift (16-P-APSK) constellation 46 illustrating 

points of the constellation of FIG. 3 in accordance with a a locus 48 of a quadrature phase-point signal stream 50 

preferred embodiment of the present invention. ( FIG * 2 ) over twelve exemplary sequential phase points 52 

in accordance with a preferred embodiment of the present 

DETAILED DESCRIPTION OF THE 60 invention. The following discussion refers to FIGS. 2 

PREFERRED EMBODIMENTS through 3. 

FIG. 1 depicts a simplified block diagram of a digital Phase mapper 44 maps symbols (i.e., binary data units) 

communications system 20 and FIG. 2 depicts a block present in interleaved signal stream 42, encoded signal 

diagram of a constrained-envelope digital communications stream 38, or input signal stream 34, into phase points 54 in 

transmitter circuit 22 in accordance with a preferred embodi- 65 phase-point constellation 46. While constellation 46 is 

ment of the present invention. The following discussion depicted in FIG. 3 as a 16-P-APSK constellation, those 

refers to FIGS. 1 and 2. skilled in the art will appreciate that the circuitry and 
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methodology of the present invention may be applied to all 
forms of constellations. The present invention is especially 
beneficial when used with constellations having rings of 
different magnitudes, i.e., amplitude and phase -shift keying 
(APSK) constellations. This is true because APSK 
constellations, requiring amplitude modulation of the signal, 
desirably use linear amplifiers to reproduce that amplitude 
modulation. 

Each phase point 54 in constellation 46 represents a 
plurality, in this example four, of symbols. The values of the 
symbols in a given phase point 54 determine the location of 
that phase point 54 within constellation 46 in a manner well 
known to those skilled in the art. 

Each quadrature phase point 54 may be thought of as 
having a vector value expressed as I,Q in the Cartesian 
coordinate system, where I is the in -phase (abscissa) value 
and Q is the quadrature (ordinate) value of the vector, or 
expressed as M,(j> in the polar coordinate system, where M 
is the magnitude and $ is the phase angle of the vector. In 
this discussion, the M,<j> designation will be used throughout, 
as the vector magnitude is the most discussed vector com- 
ponent. 

In the exemplary 16-P-APSK constellation 46 of FIG. 3, 
each phase point 54 resides upon an outer ring 56 or an inner 
ring 58. Phase -points 54 residing upon outer ring 56 are 
outer-ring or maximum-magnitude phase points 60. That is, 
outer-ring phase points 60 have a maximum magnitude 
(maximum value of M) as represented by the radius of outer 
ring 56. For purposes of discussion, the magnitudes of 
outer-ring phase points 60 are normalized to 1.00. 

Inner-ring phase points 62, i.e., those phase points 54 
residing upon inner ring 58, have a lesser magnitude as 
represented by the radius of inner ring 58. For the exemplary 
16-P-APSK constellation 46 depicted in FIG. 3, the mag- 
nitudes of inner-ring phase points 62 may desirably be 
approximately 0.63 when outer-ring phase point 60 magni- 
tudes are normalized to 1.00. 

FIG. 4 depicts a plurality of signal streams, in accordance 
with a preferred embodiment of the present invention. The 
following discussion refers to FIGS. 2 through 4. 

The output of phase mapper 44 is phase-point signal 
stream 50. Phase mapper 44 processes one phase point 54 
per unit baud interval 64. That is, phase-point signal stream 
50 consists of a series of consecutive phase-point pulses 66, 
each of which represents one phase point 54, whose leading 
edges are one unit baud interval 64 apart. Those skilled in 
the art will appreciate that other embodiments of phase-point 
signal stream 50 are equally valid, that the embodiment 
utilized is dependent upon the circuitry producing and 
processing phase-point signal stream 50, and that the use of 
other embodiments of this or any other signal stream does 
not depart from the spirit of the present invent nor the scope 
of the appended claims. 

FIGS. 3 and 4 illustrate a series of twelve exemplary 
sequential phase points 52, representative of a random data 
stream processed by transmitter circuit 22 (FIG.2). These 
twelve exemplary phase points 52 reside at temporally 
consecutive locations labeled Iq, i l7 t^, t 3 , t 4 , t 5 , t 6 , t 7 , t 8 , tg, 
t 10 , and t 1;l . These labels represent sequential integral times 
at unit baud intervals 64, i.e., integral-baud times, and 
indicate the leading-edge times of phase-point pulses 66. For 
purposes of simplification within this discussion, any occur- 
rence at time t^ shall be referred to as "occurrence tW\ For 
example, an exemplary phase point 52 occurring at time t 2 
shall be referred to as phase point t 2 , and the associated 
phase-point pulse 66 whose leading edge occurs at time t 2 
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shall be referred to as phase-point -signal pulse t2- In other 
words, at time t 2 , phase point ^ is clocked and phase-point- 
signal pulse tj begins. One unit baud interval 64 later, at time 
tj, phase point t 3 is clocked and phase-point pulse t 3 begins. 
This process continues indefinitely, with twelve exemplary 
phase points t 0 through t 21 depicted in FIG. 3 and twelve 
corresponding phase-point-signal pulses t 0 through t n 
depicted in phase-point signal stream 50 of FIG. 4. 

Table 1 below illustrates the magnitudes for phase-point- 
signal pulses to through t n . 

TABLE 1 



Phase- Point Pulse Magnitudes 


Phase-Point- Signal 




Pulse 


Magnitude 


to 


Outer-Ring 68 




Inner-Ring 70 
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Outer-Ring 68 
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Outer-Ring 68 




Inner-Ring 70 




Outer-Ring 68 




Outer-Ring 68 


h 


Outer-Ring 68 




Outer-Ring 68 




Inner-Ring 70 




Outer-Ring 68 


til 


Inner-Ring 70 



Phase point to is an outer-ring phase point 60. Phase- 
point-signal pulse to therefore has an outer- ring magnitude 
68. In like manner, phase point t 2 is an inner-ring phase point 
62 and phase-point-signal pulse t x has an inner-ring magni- 
tude 70. 

Phase-point signal stream 50 effects locus 48 through 
constellation 46. Locus 48 coincides with the location of 
each exemplary phase point ^ through t n in turn at unit baud 
intervals 64. In FIG. 3, locus 48 is depicted as effecting a 
minimum distance (straight line) path between adjacent 
exemplary phase points 52. Those skilled in the art will 
appreciate that locus 48 is so depicted solely for the sake of 
simplicity, and that in actual practice, locus 48 instantly 
jumps or snaps between exemplary phase points 52 in a 
discontinuous manner. 

FIG. 5 depicts an expanded phase-point constellation 46' 
illustrating a locus 72 of a filtered. signal stream 74 (FIG. 2) 
over twelve exemplary sequential phase points 52 in accor- 
dance with a preferred embodiment of the present invention. 
The following discussion refers to FIGS. 2 through 5. 

In the preferred embodiment, phase-point signal stream 
50 passes to the input of a pulse-spreading filter 76, pref- 
erably realized as a Nyquist-type filter, such as a Nyquist, 
root-Nyquist, raised cosine- rolloff, etc., filter. Pulse- 
spreading filter 76 filters phase -point signal stream 50 into 
filtered signal stream 74, depicted in FIG. 5. In orthogonal 
frequency division multiplex (OFDM) systems, also known 
as multitone modulation (MTM) systems, pulse-spreading 
filter 76 may be implemented using a transmultiplexer or 
equivalent circuitry. 

In accordance with Shannon's theory, well known to those 
skilled in the art, pulse-spreading filter 76 produces at least 
two (only two in the preferred embodiment) output filtered- 
signal pulses 78, i.e., complex samples of filtered signal 
stream 74, for each input phase-point pulse 66 received. This 
is demonstrated in FIG. 4 where filtered signal stream 74 
possesses two filtered-signal pulses 78 per unit baud interval 
64. In the preferred embodiment, filtered-signal pulses 78 
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consist of alternating on-time pulses 80, i.e., samples of 
filtered signal stream at integral unit baud intervals 64, and 
off- time pulses 82, i.e., samples of filtered signal stream 74 
between integral unit baud intervals. In effect, filtered signal 
stream 74 is made up of two interleaved data streams, an 
on-time signal stream 84 and an off-time signal stream 86. 

On-time signal stream 84 is substantially a version of 
phase-point signal stream 50, wherein each phase-point 
pulse 66 has been reduced in duration from one unit baud 
interval 64 to a half -unit baud interval 88 to become on-time 
pulse 80 while maintaining substantially the same relative 
leading-edge time. That is, filtered-signal pulse to has sub- 
stantially the same magnitude and substantially the same 
leading edge time as phase-point pulse to with approxi- 
mately one-half the duration. Of course, those skilled in the 
art will appreciate that signal streams 74 and 84 may be 
delayed from signal stream 50 by a delay imposed by filter 
76. 

The generation of both on-time pulses 80 and off-time 
pulses 82 by pulse-spreading filter 76 effectively populates 
expanded constellation 46' (FIG. 5) with on-time phase- 
points 90 (circles) and off-time phase points 92 (squares). 
The original phase points 54 of constellation 46 (FIG. 3), 
i.e., the phase points carrying the intelligence to be com- 
municated by transmitter circuit 22, are on-time phase points 
90 of expanded constellation 46'. 

Added to expanded constellation 46' are off- time phase 
points 92, with each off-time phase-point 92 occurring 
approximately midway in time between consecutive on-time 
phase points 90. Therefore, exemplary sequential phase 
points 52 become exemplary filtered phase points 94. Exem- 
plary filtered phase points 94 are made up of alternating 
exemplary on-time filtered phase points 96 and exemplary 
off-time filtered phase points 98, and reside at temporally 
consecutive locations labeled t 0 , t 0 5 , t 2 , t 1>5 , t^, ^ 5 , t3, t 3 5 , 

In FIG. 5, exemplary on-time filtered phase points 96 are 
located at integral-baud times (t 0 , t 1( t 2 , etc.), whereas 
exemplary off-time filtered phase points 98 are located at 
fractional-baud (non-integral-baud) times (io, S9 t ltS , t 2 5 , etc.). 

The generation of off-time phase points 92 approximately 
midway in time between consecutive on-time phase points 
90 causes filtered signal locus 72 to effect excursions having 
local peak magnitudes 99 greater than outer-ring magnitude 
68. Such excursions occur because the immediate position 
of locus 72 at any given instant in time is not only a result 
of those phase points 54 proximate that position, but of a 
plurality of phase points 54 both preceding and following 
that instant in time. That is, in the preferred embodiment, the 
determination of the position of locus 72 at time t 2 s (i.e., 
coincident with off-time phase point ^ 5 ) is determined not 
only by the positions of phase points t^ and l,, but by the 
positions of numerous phase points 54 preceding phase point 

5 (i.e., phase points t^ t 1>5 , X 19 {q 5 , etc.) and the positions 
of numerous phase points 54 following phase point ^ 5 (i.e., 
phase points t^, t 3 5 , t 4 , t 4 5 , etc.) 

This phenomenon is illustrated in FIG. 6, which depicts a 
pair of Nyquist-type datum bursts 100 in accordance with a 
preferred embodiment of the present invention. The follow- 
ing discussion refers to FIGS. 2, 4, 5, and 6. 

In the preferred embodiment, pulse -spreading filter 76 is 
realized as a Nyquist-type filter. Therefore, when a single 
phase-point pulse 66 is filtered by pulse-spreading filter 76, 
that single pulse 66 is transformed into a Nyquist-type 
datum burst 100 extending over a plurality of unit baud 
intervals 64. It is a property of Nyquist-type filters that 
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datum burst 100 attains a datum-burst peak value 102 (i.e., 
a local peak magnitude) at the primary sampling time of the 
specific phase-point pulse 66 (i.e., at time t 2 for phase-point 
pulse Q, and attains a zero datum-burst value 104 (i.e., is 

5 equal to zero) at integral unit baud intervals 64 preceding 
and following peak datum-burst value 102 (i.e., at times . . . , 
t_ a , t 0 , tj, and t 3 , t 4 , t 5 , . . . , for phase point pulse tj. In this 
manner, the energy of each pulse 78 is spread over a plurality 
of baud intervals 64 preceding and following the clocking 

1Q instant (time tj). 

FIG. 6 illustrates Nyquist-type datum bursts 100 for 
phase-point pulses ^ and t3, with datum burst t 2 depicted as 
a solid line and datum burst t 3 depicted as a dashed line. As 
an example, it may be seen from FIG. 6 that at time t^ the 

35 value of datum burst 12 is peak datum-burst value 102. At 
every other time separated from time t 2 by an integral 
number of unit baud intervals 64, the value of datum burst 
is zero. An analogous condition occurs for datum burst t 3 . 
The value of locus 72 is, at each moment in time, the sum 

2Q of all datum bursts 100 at that moment. In the simplified 
two-datum-burst example of FIG. 6, locus 72, depicted by a 
dotted line, is the sum of datum burst ^ and datum burst Lj. 
Since datum bursts ^ and t 3 are zero at each integral time t^ 
except times ^ and t 3 , the value of locus 72 is also zero 

25 except at times t 2 and t 3 , where it assumes the peak values 
of datum bursts t 2 and t 3 , respectively. 

The value of locus 72 at any instant in time between 
integral-baud times is the sum of the values of all datum 
bursts 100 at that instant. For example, in FIG. 6 where only 

3 0 two datum bursts 100 are considered, locus 72 has a value 
at time ^ 5 that is the sum of the values of datum bursts t 2 
and t 3 at time t2. 5 . Since datum bursts ^ and t, both have 
significant positive values at time ^ 5 , locus 72 has a value 
significantly greater than the maximum values of either 

35 datum burst t 2 or datum burst t 3 . 

Since locus 72 describes the sum of all datum bursts 100, 
locus 72 is a function of the shape of the curves (FIG.6) 
describing those datum bursts 100. That is, locus 72 is a 
function of a filtered-signal peak magnitude component of a 

40 filtered-signal complex digital value at any given point. The 
shape of the datum-burst curve is a function of the excess 
bandwidth factor, a, a design property of pulse-spreading 
filter 76. The smaller the value of a, the more locus 72 may 
rise above the peak datum burst values 102 of adjacent 

45 datum bursts 100. Typical designs of pulse-spreading filters 
76 use a values of 0.15 to 0.5. For like-valued adjacent 
phase points 54 and an a value of 0.2, a maximum excursion 
magnitude 105 (i.e., the potential local peak magnitude 99 of 
locus 72) is approximately 1.8 times the value of the 

50 maximum phase -point magnitude. That is, the magnitude of 
the constrained envelope is approximately 1.8 times that of 
the unconstrained envelope. In the preferred embodiment 
depicted in FIGS. 3, 4, and 6, on-time phase points ^ and t, 
are both outer-ring phase points 60 having a normalized 

55 outer-ring magnitude 68 of 1.00. Therefore, off-time phase 
point ^ 5 may have a normalized maximum excursion mag- 
nitude 105 of 1.8. This implies that transmitter circuit 22, to 
faithfully transmit phase point t 25 without excessive 
distortion, and without the benefit of the present invention, 

60 would require an output power of 3.24 (L8 2 ) times the 
power required to transmit phase point t 2 or 1$, which are 
representative of the highest magnitude intelligence- 
carrying phase points 54. This represents an inefficient use 
of available power. 

65 The following discussion refers to FIGS. 2, 4, and 5. 
Off-time signal stream 86, a portion of filtered signal 
stream 74, passes from an output of pulse-spreading filter 76 
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to an input of an off-time constrained -envelope generator 
106. It is the task of off-time constrained-envelope generator 
106 to produce an off-time constrained-bandwidth error 
signal stream 108 from off-time signal stream 86. A complex 
summing or combining circuit 110 combines off-time 
constrained-bandwidth error signal stream 108 with a 
delayed version of filtered signal stream 74 (discussed 
below) to produce a constrained-envelope signal stream 112. 
Constrained -envelope signal stream 112 is effectively fil- 
tered signal stream 74 with compensation for excursions of 
locus 72 with magnitudes greater than outer-ring magnitude 
68. 

A quadrature threshold generator 118 generates a quadra- 
ture threshold signal 120. In the preferred embodiment, 
threshold signal 120 is a steady-state, constant signal having 
a value approximately equal to outer-ring magnitude 68. 
Threshold signal 120 is used to establish a reference with 
which off-time signal stream 86 is compared. Those skilled 
in the art will appreciate that threshold signal 120 may 
assume many forms and values in keeping with the meth- 
odology and circuitry incorporated in the comparison. The 
use of other forms and/or other values does not depart from 
the spirit of the present invention nor from the scope of the 
appended claims. 

Threshold signal 120 and off-time signal stream 86 are 
combined in an off-time complex summing or combining 
circuit 122 to produce an off-time difference signal stream 
124. Off-time difference signal stream 124 is made up of a 
series of off-time difference pulses 126 whose values are the 
difference between the values of equivalent off-time pulses 
82 and the value of threshold signal 120. Since any given 
off-time pulse 82 may have a value greater than, equal to, or 
less than the value of threshold signal 120, off-time differ- 
ence signal stream 124 would normally be made up of a 
combination of off-time difference pulses 126 having 
positive, zero, and negative values. 

Off-time difference signal stream 124 is passed to the 
input of an off-time discriminator 128 to produce an off-time 
error signal stream 130. In the preferred embodiment, off- 
time error signal stream 130 is a variation of off-time 
difference signal stream 124 in which all off-time difference 
pulses 126 having positive values are passed unchanged as 
off-time error pulses 132 while all other off-time difference 
pulses 126 are passed as zero-value pulses (i.e., eliminated). 
In other words, off-time error signal stream 130 is formed 
from pulses, the timing of which coincide with excursions of 
locus 72 beyond outer-ring magnitude 68 and the magni- 
tudes of which correspond to the degree to which locus 72 
passes beyond outer-ring magnitude 68. 

Off-time error signal stream 130 is then passed to the 
input of an off-time pulse-spreading filter 134. Off-time 
pulse -spreading filter 134 is substantially identical to first 
pulse-spreading filter 76. That is, in the preferred 
embodiment, both pulse spreading filters 76 and 134 are 
realized as Nyquist-type filters with substantially identical 
transfer characteristics. Off-time pulse-spreading filter 134 
produces off-time constrained-bandwidth error signal stream 
108 and completes the action of off- time constrained- 
envelope generator 106. 

Within off-time constrained-envelope generator 106, off- 
time pulse-spreading filter 134 receives one off-time error 
pulse 132 from off-time discriminator 128 per unit baud 
interval 64. Off-time pulse -spreading filter 134 then trans- 
forms each off-time error pulse 132 into a Nyquist-type error 
burst (not shown) extending over a plurality of unit baud 
intervals. Since off-time pulse -spreading filter 134 is a 
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Nyquist-type filter, each error burst attains an error-burst 
peak value (not shown) at the primary sampling time of the 
specific off-time error pulse 132 (i.e., at time ^ 5 for error 
pulse 5 ), and attains a zero error-burst value (not shown) 

5 at integral unit baud intervals 64 preceding and following 
the peak error-burst value (i.e., at times . . . , \_ x 5 , Iq s , t 1-5 , 
and tg 5 , t 4 5 , t 5 5 , . . . , for error pulse tj 5 ). In this manner, 
the energy of each off-time constrained -envelope error pulse 
136 is spread over a plurality of baud intervals 64 preceding 

1Q and following the clocking instant (time ^ 5 ). This results in 
the conversion of off-time error signal stream 130 into 
off-time constrained-bandwidth error signal stream 108. 
Off- time constrained -bandwidth error signal stream 108 is 
made up of off-time constrained-envelope error pulses 136. 
This operation is essentially the same as the operation of 

15 pulse-spreading filter 76 in the conversion of phase-point 
signal stream 50 into filtered signal stream 74 described 
hereinabove. 

Since off-time constrained-envelope error pulses 136 are 

2Q derived from off-time pulses 82, the error-burst peak and 
zero values occur approximately midway between integral 
baud times, i.e., at baud times t^, t a 5 , t 25 , etc., hence 
between datum-burst peak and zero values 102 and 104 of 
filtered signal stream 74. 

25 The production of off-time constrained-bandwidth error 
signal stream 108 completes the operation of off-time con- 
strained envelope generator 106. 

Filtered signal stream 74 is also passed to the input of a 
delay element 138. Delay element 138 produces delayed 

30 signal stream 140, which is effectively filtered signal stream 
74 delayed sufficiently to compensate for the propagation 
and other delays encountered in off-time constrained- 
envelope generator 106, and particularly in off-time pulse- 
spreading filter 134. In other words, delayed signal stream 

35 140 is filtered signal stream 74 brought into synchronization 
with off-time constrained-bandwidth error signal stream 
108. 

Combining circuit 110 combines filtered signal stream 74, 
in the form of delayed signal stream 140, and off-time 

40 constrained-bandwidth error signal stream 108 to reduce 
peak magnitude components of filtered signal stream 74. A 
resultant constrained-envelope signal stream 112 is made up 
of a series of digital pulses 142 whose values are the 
difference between the values of corresponding filtered- 

45 signal pulses 78 and off-time constrained-envelope error 
pulses 136. The result is a series of digital pulses 142 whose 
values do not appreciably exceed outer-ring magnitude 68 of 
expanded constellation 46'. 

In some embodiments of the present invention, certain of 

50 outer-ring phase points 60 may have magnitudes greater 
than outer-ring magnitude 68, i.e., may be located beyond 
outer ring 56. This condition may occur as a result of 
pulse -spreading filter 76 executing certain Nyquist-type 
functions well known to those skilled in the art. In such an 

55 embodiment, transmitter circuit 22 contains an on-time 
constrained envelope generator 106' in addition to off- time 
constrained-envelope generator 106 discussed above. 

On-time signal stream 84, also a portion of filtered signal 
stream 74, passes from an output of pulse-spreading filter 76 

60 to an input of on-time constrained-envelope generator 106'. 
It is the task of on-time constrained-envelope generator 106' 
to produce an on-time constrained-bandwidth error signal 
stream 108' from on-time signal stream 84. Combining 
circuit 110 combines both off-time and on-time constrained- 

65 bandwidth error signal streams 108 and 108' with the 
delayed version of filtered signal stream 74 (discussed 
below) to produce constrained-envelope signal stream 112. 
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On-time constrained-envelope generator 106' operates in on-time phase points 90. The centers of phase-point prob- 
a manner analogous with the operation of off- time abilities 144 occupy the same normalized locations within 
constrained-envelope generator 106. Threshold signal 120 noise-influenced constellation 46" as do on-time phase 
and on-time signal stream 84 are combined in an on-time points 90 within expanded constellation 46'. 
complex summing or combining circuit 122' to produce an 5 The positional aberrations of transmitted phase points 145 
on-time difference signal stream 124*. On-time difference relative to the corresponding on-time phase points 90 rep- 
signal stream 124' is passed to the input of an on-time re sent a degree of positional error. This positional error 
discriminator 128' to produce an on-time error signal stream degrades the bit error rate and effects a detriment to trans- 
130'. On-time error signal stream 130' is then passed to the mission. The absence of off-time phase points 92 with a 
input of an on-time pulse-spreading filter 134', which pro- 10 magnitude significantly greater than outer-ring magnitude 
duces on-time constrained bandwidth error signal stream * 8 ( nG * *> 1D constramed^nvelope signal stream 112, 

108'. Like off-time pulse-spreading filter 134, on-time pulse- 7^ allo ; s an mcre ^ 10 P owcr ou ^ ut for a ^ n 

i • cu --J, . u* **Tt *■ i. a * i bandwidth and power amplifier that more than compensates 

spreading niter 134 , is substantially identical to first pulse- c +u cJ *.« j i_ • * lie a 

j- cu lor the position error of transmitted phase points 145. A net 

spreading filter 76. u 

r & improvement m performance results. 

Since on-time constrained-envelope error pulses (not ^ Referring back to FIG. 2, the output of combining circuit 

shown) are derived from on-time pulses 80, the error-burst 110 CO nstrained-envelope signal stream 112, is passed to an 

peak and zero values occur at integral baud times, i.e., at input of a subst antially linear amplifier 146. Substantially 

baud times t a , ^ t 3 , etc., hence between datum-burst peak Unear amplifier 146 pro duces RF broadcast signal 26, which 

and zero values 102 and 104 of filtered signal stream 74. ^ then broadcast ^ transmitter antenna 24. In the preferred 

Combining circuit 110 combines filtered signal stream 74, 20 embodiment, substantially Unear amplifier 146 is made up of 

in the form of delayed signal stream 140, with both off- time a digital linearizer 148, a digital-to-analog converter 150, 

and on-time constrained-bandwidth error signal stream 108 and a radio-frequency (RF) amplifying circuit 152. Those 

and 108' to reduce peak magnitude components of filtered skilled in the art will appreciate that substantially linear 

signal stream 74. ^ amplifier 146 may be realized in any of a plurality of 

A side effect of this methodology is that locus 72 at different embodiments other than that described here, and 

integral unit baud intervals 64 adds a signal-dependent, that utilization of any of these different embodiment does not 

baud -limited noise factor to the positions of phase points 54 depart from the intent of the present invention nor the scope 

in constellation 46 (FIG. 3). This results in transmitter circuit of the appended claims. 

22 transmitting a "noise-influenced" phase-point constella- 3Q Within substantially linear amplifier 146, digital linearizer 
tion 46". In FIG. 7, noise-influenced constellation 46" is 148 alters constrained -envelope signal stream 144 into a 
depicted illustrating constrained-envelope phase -point prob- pre-distorted digital signal stream 154. Pre-distorted digital 
abilities 144 of phase points 54 in accordance with a signal stream 154 is made non-linear in just the right manner 
preferred embodiment of the present invention. The follow- to compensate for non-linearities within digital-to-analog 
ing discussion refers to FIGS. 2, 3, 5 and 7. 35 converter 150 and RF amplifying circuit 152, hence linear- 
Phase -point probabilities 144 reside in noise-influenced izing substantially linear amplifier 146. 
constellation 46" exactly as phase points 54 reside in con- Digital-to- analog converter 150 then converts pre- 
ste Uation 46, i.e., in the same configuration with centers at distorted digital signal stream 154 into an analog baseband 
the same locations. The actual location of a given transmit- signal 156. Analog baseband signal 156 is then amplified by 
ted phase point 145 within a given phase-point probability 40 RF amplifying circuit 152 into RF broadcast signal 26 and 

144 is a function of a plurality of variable conditions and, transmitted via transmitter antenna 24. 

although somewhat correlated, except in certain specialized In summary, the present invention teaches a methodology 

cases, cannot readily be predicted. In effect, for a given and circuitry by which a transmitter circuit utilizing 

phase point 54, the resultant transmitted phase point 145 Nyquist-type filtration may produce a constrained envelope 

may be located anywhere within phase-point probability 45 having a magnitude at or near the approximate uncon- 

144, i.e., within an indeterminate area having a center strained envelope magnitude of the desired constellation, 

coincident with the location of the original phase point 54. This enables the transmitter output amplifier to be biased so 

The probability of transmitted phase point 145 being located that the maximum unconstrained envelope magnitude is at 

at any specific position within that indeterminate area varies or near the top of the amplifier's linear region without 

as an inverse function of the distance of that specific position 50 incurring clipping of the constrained envelope transmis- 

from the location of the original phase point 54. sions. This in turn produces a more efficient output amplifier 

For any given phase point 54, the transmitted phase point and effects an increase in the power output of a given output 

145 may be said to be proximate its idealized position within amplifier. Conversely, a lower power amplifier may be used 
noise -influenced constellation 46". That is, a locus (not to provide the same output power that was previously 
shown) of constrained-envelope signal stream 112 passes 55 output. This effects a significant savings in output amplifier 
proximate the idealized positions of exemplary phase points cost. 

to» t 1? t2, etc., at the clocking instants in time. Although the preferred embodiments of the invention 

The original phase points 54 of constellation 46, as have been illustrated and described in detail, it will be 

produced by phase mapper 44, are on-time phase points 90 readily apparent to those skilled in the art that various 

(circles) of expanded constellation 46'. It is these on-time 60 modifications may be made therein without departing from 

phase points 90 that carry the intelligence of RF broadcast the spirit of the invention or from the scope of the appended 

signal 26 as ultimately transmitted. Off-time phase points 92 claims, 

(squares) are by-products of pulse-spreading filter 76, What is claimed is: 

required to constrain spectral re growth, and carry no intel- 1. A constrained-envelope digital communications trans- 

ligence. Phase-point probabilities 144 of noise-influenced 65 mitter circuit comprising: 

constellation 46" represent the resultant areas of probable a pulse -spreading filter configured to receive a quadrature 

locations of transmitted phase points 145 as derived from phase -point signal stream of digitized quadrature phase 



12/30/04, EAST Version: 2.0.1.4 



6,104,761 



13 



14 



points and produce a filtered signal stream, said filtered 
signal stream exhibiting energy corresponding to each 
phase point spread throughout a plurality of unit baud 
intervals; 

a constrained-envelope generator coupled to said pulse- 
spreading filter and configured to produce a 
constrained-bandwidth error signal stream; 

a combining circuit coupled to said pulse-spreading filter 
and to said constrained-envelope generator, said com- 
bining circuit configured to combine said filtered signal 
stream and said constrained-bandwidth error signal 
stream to produce a constrained-envelope signal 
stream; and 

a substantially linear amplifier having an input coupled to 
said combining circuit. 

2. A digital communications transmitter circuit as claimed 
in claim 1 wherein said pulse-spreading filter is a Nyquist- 
type filter. 

3. A digital communications transmitter circuit as claimed 
in claim 1 wherein said combining circuit is configured to 
combine said filtered signal stream and said constrained- 
bandwidth error signal stream to reduce a peak magnitude 
component of said filtered signal stream. 

4. A digital communications transmitter circuit as claimed 
in claim 3 wherein said combining circuit is a complex 
summing circuit. 

5. A digital communications transmitter circuit as claimed 
in claim 1 wherein: 

said pulse-spreading filter is a first pulse-spreading filter; 

said transmitter circuit additionally comprises a delay 
element coupled between said first pulse-spreading 
filter and said combining circuit; and 

said constrained-envelope generator comprises a second 
pulse -spreading filter coupled to said combining cir- 
cuit. 

6. A digital communications transmitter circuit as claimed 
in claim 5 wherein: 

said first pulse -spreading filter is configured so that each 
phase point is transformed into a Nyquist-type datum 40 
burst extending over a plurality of unit baud intervals, 
having a datum -burst peak value occurring in one of 
said plurality of unit baud intervals and datum-burst 
zero values occurring substantially at integral unit baud 
intervals away from said datum-burst peak value, so 45 
that said filtered signal stream in each unit baud interval 
substantially equals the sum of said Nyquist-type 
datum bursts from a plurality of phase points; and 

said constrained-envelope generator is configured so that 
said second pulse -spreading filter receives error pulses, 50 
transforms each error pulse into a Nyquist-type error 
burst extending over a plurality of unit baud intervals, 
having an error-burst peak value occurring in one of 
said plurality of unit baud intervals and error-burst zero 
values occurring substantially at integral unit baud 
intervals away from said error-burst peak value, so that 
said constrained -bandwidth error signal stream in each 
unit baud interval substantially equals the sum of said 
Nyquist-type error bursts from a plurality of error 
pulses. 

7. A digital communications transmitter circuit as claimed 
in claim 6 wherein said constrained -envelope generator is 
configured so that said Nyquist-type error bursts exhibit said 
error-burst peak values and said error-burst zero values at 
instances in time when said Nyquist-type datum bursts 
exhibit neither said datum-burst peak values nor said datum- 
burst zero values. 
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8. A digital communications transmitter circuit as claimed 
in claim 7 wherein said constrained-envelope generator is 
configured so that said error-burst peak values and said 
error-burst zero values occur approximately midway 
between said datum-burst peak values and said datum-burst 
zero values. 

9. A digital communications transmitter circuit as claimed 
in claim 5 wherein said first and second pulse-spreading 
filters exhibit substantially equivalent transfer characteris- 
tics. 

10. A digital communications transmitter circuit as 
claimed in claim 5 wherein: 

said first pulse -spreading filter receives one quadrature 
phase point per unit baud interval and produces two 
complex samples of said filtered signal stream per unit 
baud interval; 

said constrained-envelope generator evaluates one of said 
two complex samples of said filtered signal stream 
produced by said first pulse-spreading filter per unit 
baud interval; and 

said second pulse-spreading filter receives one error pulse 
per unit baud interval and produces two complex 
samples of said constrained-envelope error-signal 
stream per unit baud interval. 

11. A digital communications transmitter circuit as 
claimed in claim 1 wherein: 

said filtered signal stream is a stream of complex digital 
values, with each of said complex digital values exhib- 
iting a peak magnitude component; and 

said constrained-envelope generator is configured to 
determine when ones of said peak magnitude compo- 
nents exceed a threshold value. 

12. A digital communications transmitter circuit as 
claimed in claim 11 wherein: 

said transmitter circuit additionally comprises a phase 
mapper coupled to said pulse-spreading filter and con- 
figured to select said digitized quadrature phase points 
from a phase-point constellation, said phase-point con- 
stellation having a maximum-magnitude phase point; 
and 

said threshold value is a magnitude value approximately 
equal to a magnitude of said maximum-magnitude 
phase point. 

13. A digital communications transmitter circuit as 
claimed in claim 1 additionally comprising an interleaver 
coupled to said phase mapper. 

14. A digital communications transmitter circuit as 
claimed in claim 1 wherein: 

said constrained-envelope generator is an off-time 
constrained-envelope generator; 

said constrained-bandwidth error signal stream is an off- 
time constrained-bandwidth error signal stream; 

said transmitter circuit additionally comprises an on-time 
constrained-envelope generator coupled to said pulse - 
spreading filter and configured to produce an on-time 
constrained-bandwidth error signal stream; and 

said combining circuit is coupled to said pulse-spreading 
filter, to said off-time constrained-envelope generator, 
and to said on- time constrained-envelope generator, 
and said combining circuit is configured to combine 
said filtered signal stream, said off-time constrained- 
bandwidth error signal stream, and said on-time 
constrained-bandwidth error signal stream to produce 
said constrained-envelope signal stream. 

15. A digital communications transmitter circuit as 
claimed in claim 1 wherein said substantially linear ampli- 
fier comprises; 



12/30/04, EAST version: 2.0.1.4 



6,104,761 



15 



16 



15 



20 



a digital linearizer configured to pre-distort said 
constrained-envelope signal stream into a pre-dislorted 
digital signal stream; 

a digital-to-analog converter coupled to said digital lin- 
earizer and configured to produce an analog baseband 5 
signal from said pre-distorted digital signal stream; and 

a radio- frequency amplifying circuit configured to gener- 
ate a radio-frequency broadcast signal from said analog 
baseband signal. 

16. In a digital communications system, a method for the 10 
transmission of a constrained-envelope communications 
signal, said transmission method comprising the steps of: 

filtering a quadrature phase-point signal stream to produce 
a filtered signal stream, said filtering step spreading 
energy from each phase point in said filtered signal 
stream over a plurality of unit baud intervals; 

generating a constrained-bandwidth error signal stream 
from said filtered signal stream and a threshold signal; 

combining said filtered signal stream and said 
constrained-bandwidth error signal stream to produce a 
constrained-envelope signal stream; 

linearly amplifying said constrained-envelope signal 
stream to produce said constrained-envelope commu- 
nications signal; and 

transmitting said constrained -envelope communications 
signal. 25 

17. A transmission method as claimed in claim 16 wherein 
said combining step comprises the step of reducing a peak 
magnitude component of said filtered signal stream. 

18. A transmission method as claimed in claim 16 
wherein: 30 

said generating step comprises the step of filtering an 
error signal stream having one error pulse per unit baud 
interval to produce said constrained-bandwidth error 
signal stream, said filtering step spreading energy from 
each error pulse in said error signal stream over a 35 
plurality of unit baud intervals; 

said transmission method additionally comprises the step 
of delaying said filtered signal stream to produce a 
delayed signal stream; and 

said combining step combines said delayed signal stream 40 
and said constrained-bandwidth error signal stream to 
produce said constrained-envelope signal stream. 

19. A transmission method as claimed in claim 16 
wherein: 

said filtering step comprises the step of receiving one 45 
quadrature phase point per unit baud interval; 

said filtering step additionally comprises the step of 
producing two complex samples of said filtered signal 
stream per unit baud interval; 5Q 

said generating step comprises the step of evaluating one 
of said two complex samples of said filtered signal 
stream per unit baud interval to produce an error signal 
stream having one error pulse per unit baud interval; 
and 55 

said generating step additionally comprises the step of 
filtering said error signal stream to produce said 
constrained-bandwidth error signal stream having two 
complex samples of said constrained -band width error 
signal stream per unit baud interval. 60 

20. A transmission method as claimed in claim 19 wherein 
said generating step additionally comprises the steps of: 

providing said threshold signal; and 

determining when ones of peak magnitude components of 
a stream of complex digital values of said filtered signal 65 
stream exceed a threshold value of said threshold 
signal 



21. A transmission method as claimed in claim 16 
wherein: 

said filtered signal stream includes two or more complex 
digital values per unit baud interval, said complex 
digital values in said filtered signal stream exhibiting 
local peak magnitudes; and 

said generating step is configured so that said constrained - 
bandwidth error signal stream includes two or more 
complex values per unit baud interval, said complex 
values in said constrained-bandwidth error signal 
stream being responsive to said local peak magnitudes 
of said filtered signal stream so as to spread energy 
from selected ones of said local peak magnitudes over 
a plurality of unit baud intervals of said constrained- 
bandwidth error signal stream. 

22. A transmission method as claimed in claim 16 wherein 
said transmitting step continuously transmits said 
constrained-envelope communications signal. 

23. A constrained-envelope digital -communications trans- 
mitter circuit comprising: 

a binary data source configured to provide a binary input 
signal stream; 

a phase mapper coupled to said binary data source and 
configured to produce a quadrature phase-point signal 
stream, wherein said phase-point signal stream has a 
predetermined number of symbols per unit baud 
interval, said predetermined number of symbols defin- 
ing a phase point in a phase-point constellation; 

a Nyquist-type filter coupled to said phase mapper and 
configured to produce a filtered signal stream; 

a constrained-envelope generator coupled to said 
Nyquist-type filter and configured to produce a 
constrained-bandwidth error signal stream; 

a delay element coupled to said Nyquist-type filter and 
configured to produce a delayed signal stream synchro- 
nized with said constrained-bandwidth error signal 
stream; 

a complex summing circuit coupled to said delay clement 
and said constrained-envelope generator and config- 
ured to produce a constrained-envelope signal stream; 
and 

a substantially linear amplifier coupled to said complex 
summing circuit and configured to produce a radio- 
frequency broadcast signal. 

24. A digital-communications transmitter circuit as 
claimed in claim 23 wherein said Nyquist-type filter is a first 
Nyquist-type filter, said filtered signal stream includes a first 
filtered-signal data stream and a second filtered-signal data 
stream, and said complex summing circuit is a first complex 
summing circuit, wherein said transmitter circuit addition- 
ally comprises a quadrature threshold generator configured 
to-provide a threshold signal, said threshold signal having a 
threshold value, and wherein said constrained-envelope gen- 
erator comprises: 

a complex summing circuit coupled to said first Nyquist- 
type filter and said quadrature threshold generator and 
configured to produce a difference 1 signal stream, 
wherein said difference signal stream is a stream of 
difference pulses having difference-pulse values of a 
first polarity and difference-pulse values of a second 
polarity; 

a discriminator coupled to said complex summing circuit 
and configured to produce an error signal stream from 
said difference signal stream, wherein said error signal 
stream is a stream of error pulses substantially coinci- 
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dent with said difference pulses of said difference 
signal stream, and wherein, when ones of said differ- 
ence pulses have said first-polarity difference-pulse 
values, said coincident error pulses have error-pulse 
values substantially equal to said first-polarity 5 
difference-pulse values, and when ones of said differ- 
ence pulses have said second-polarity difference-pulse 
values, said coincident error pulses have error-pulse 
values substantially equal to zero; and 
a second Nyquist-type filter coupled to said discriminator 10 
and configured to produce said constrained-bandwidth 
error signal stream. 

25. A digital-communications transmitter circuit as 
claimed in claim 24 wherein said transmitter circuit addi- 
tionally comprises: 15 

a convolutional encoder coupled to said binary data 
source and configured to produce an encoded signal 
stream; and 

an interleaver coupled to said convolutional encoder and 2Q 
configured to produce an interleaved signal stream by 
temporally decorrelating said encoded signal stream. 

26. A digital-communications transmitter circuit as 
claimed in claim 24 wherein: 

said filtered signal stream is a quadrature signal stream 
having a locus that passes proximate one of said phase 



points of said phase-point constellation at integral unit 
baud intervals; 
said first filtered-signal data stream comprises on-time 
samples of said filtered signal stream, each of said 
on-time samples occurring substantially coincident ally 
with said passage of said filtered signal locus proximate 
one of said phase points of said phase-point constella- 
tion; and 

said second filtered-signal data stream comprises off-time 
samples of said filtered signal stream wherein each of 
said off-time samples occurs between adjacent ones of 
said on-time samples. 

27. A digital-communications transmitter circuit as 
claimed in claim 26 wherein each of said off-time samples 
occurs substantially midway between adjacent ones of said 
on-time samples. 

28. A digital-communications transmitter circuit as 
claimed in claim 23 additionally comprising an interleaver 
coupled to said binary data source and configured to provide 
an interleaved signal stream. 

29. A digital-communications transmitter circuit as 
claimed in claim 23 wherein said constellation is an ampli- 
tude and phase shift keying constellation. 
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